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THE  MEASUREMENT  OF  NON-LINEAR  FORCES  AND  MOMENTS 
BY  MEANS  OF  FREE  FLIGHT  TESTS 

ABSTRACT 

It  has  been  observed,  that  the  behavior  of  missiles  either 
moving  under  the  influence  of  non-linear  forces  and  moments  or  flying 
at  large  angles  of  yaw  is  frequently  well  described  by  curves  of  the 
same  form  as  those  generated  by  linear  force  systems  and  small  angles 
of  yaw..  With  this  in  mind  an  "equivalent"  linear  solution  to  the 
actual  equations  of  yawing  motion  is  obtained. 

This  equivalent  linear  solution  has  been  used  in  the  analysis 
of  a  wide  variety  of  programs  fired  on  BRL's  Spark  Ranges  and  considerable 
success  has  been  experienced.  Excellent  internal  consistency  has 
been  observed  in  measuring  non-linear  normal  and  Magnus  forces  and 
their  moments  and,  in  all  cases  where  wind  tunnel  results  were  available, 
they  were  in  good  agreement  with  range  results. 

The  application  of  this  technique  to  the  equally  important  problem 
of  predicting  yawing  motion  is  described. 
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TABLE  OF  SYMBOLS* 


axial  moment  of  inertia 


2  2 
K10  +K20 


coefficients  in  the  drag  equation  (Equation  (22)) 


transverse  moment  of  inertia 


K10  K20 
diameter 


->  ->  -> 

ei>  e2>  e3 


unit  vectors  along  axes  of  fixed^plane  co-ordinates 


(Fr,  Fg,  F^)  aerodynamic  force  vector 
drag  force 


lift  force 


^(p)  * 


acceleration  due  to  gravity 
iJL."  JD  +  k2  ^JH  “  ^JMA^ 


P  <1 


8  dr  dq 


— K 


k  =  0 


Ki  2k  5  '  1  =  A>  D'  DA>  F>  H>  L'  M>  m>  N'  WAj  S>  T 

8 


(Expansion  as  a  function  of  8  of  aerodynamic  coefficients 
defined  in  Equations  (5  -  10)) 


Only  those  symbcjls  which  appear  in  the  body  of  this  report  are 
listed  here.  Symbols  which  are  introduced  in  the  appendices  appear 
close  to  their  definitions. 
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-  up 

Ki0  e  ;  i  =  1,  2  amplitude  of  i-th  frequency 

i  =  1,2  mid-range  amplitude  of  i-th  frequency  (p  =  0). 

axial  spin  deceleration  coefficients  < 

drag  coefficient 

axial  drag  coefficient 

Magnus  force  coefficient 

moment  coefficient  due  to  cross  angular  velocity 
lift  force  coefficient  due  to  yaw 

overturning  (or  righting)  moment  coefficient  due  to  yaw 

moment  coefficient  due  to  cross  acceleration 

normal  force  coefficient  due  to  yaw 

normal  force  coefficient  due  to  cross  acceleration 


k  =  0 


k  =  0 


t 


m 


(Mx,  M2,  M^)  aerodynamic  moment  vector 


mass 


n2J  )  direction  cosines  of  missile’s  axis  with  respect  to 
range  co-ordinates 

r*  u 

P  =  I  dt  arclength  along  the  trajectory  in  calibers 

-*  o 


Q-. 


constants  in  the  swerve  equations  (Equations  (82  -83)) 
change  in  center  of  mass  location  measured  in  calibers  (Eqs. 

(99  -loo)) 


T 

t 

t 

u 

u 


^JL  "  kl  JT)  =  i-x'  T; 


k  =  0 


*  2k 
2k  6 


time 


(u^,  u^,  u^)  velocity  vector 
|  "u  |  magnitude  of  the  velocity 


magnitude  of  the  velocity  at  mid-range 
(xy  x^>  Xj)  vector  in  the  range  co-ordinate  system 
(y^>  V2>  3^)  vector  in  the  fixed-plane  co-ordinate  system 

exponential  damping  coefficient  of  i-frequency 
=  \J W  magnitude  of  the  sine  of  the  total  yaw  angle 


a. 


K  +  IK  "  + 
10  +|A 20 


2  .  K10  K  '  -  K20  4’ 


^  -  0a’ 


5 


& 


el 

2 

e2 


K102  + 


K  ^  + 
20 


2K, 


20 


2K. 
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Effective  squared 
yaws 


£ 


I 


maximum  value  of  6 


orientation  angle  of  plane  of  yaw 


Ug  +  i^z 

\2  +  iX,^  =  - — - -  complex  yaw 

(au  +  i<s,)d 

^2  +  i|i3  =  - - - - —  dimensionless  complex  cross  angular 


u 


velocity 
u)^d 


u 


dimensionless  missile  spin 


A 

B  V 
fid 

— -  dimensionless  co-ordinate  system  spin 

u 

average  value  of  v 
air  density 

0iO  +  0i'  p  phase  angle  of  i-frequency 

0r'  -  -  02 ' 

1*10  ‘  ®  -  1*20 


arccos  ru 


fij)  angular  velocity  vector  of  co-ordinate  system 
cd2^  o)j)  angular  velocity  vector  of  missile 


—  (  ) 
dt  k  ‘ 


dp 


(  ) 


circumflex  superscript  with  the  exception  of  v  denotes 
quantities  in  the  fixed-plane  co-ordinate  system 


tilde  superscript  with  the  exception  of  6  denotes  quantities 
appearing  in  the  solution  of  the  linearized  yaw  equation 


asterisk  denotes  the  modification  of  quantities  involving 
aerodynamic  coefficients  through  the  consideration  of  the 
cosine  of  the  yaw  angle,  &. 


INTRODUCTION  i 


1. 

The  free  flight  spark  range  technique  measures  the  aerodynamic 
forces  and  moments  acting  on  a  missile  hy  means  of  very  accurate 
observations  of  its  motion  in  flight.  This  process  requires  a 
knowledge  of  the  functional  dependence  of  these  forces  and  moments 
on  the  dynamic  variables  of  the  motion  in  order  that  the  solution 
curves  to  the  differential  equations  of  motion  may  be  obtained. 

These  solution  curves  are  fitted  to  the  motion  and  the  forces 
and  moments  calculated  from  the  parameters  of  the  fit. 

This  need  for  solutions  in  closed  form  has  traditionally 
limited  range  tests  to  motions  which  are  described  by  linear  equations. 
Since  non-linear  terms  arise  from  both  the  size  of  the  motion  (non¬ 
linear  geometry  terms)  and  the  presence  of  second  order  or  higher 
terms  in  the  aerodynamic  force  expansion  (non-linear  force  terms), 
this  means  that  the  range  technique  is  restricted  to  configurations, 
possessing  linear  force  systems  and  flying  at  small  angles  of  yaw.. 

Strangely  enough,  a  number  of  models  have  been  fired  in  the  BRL 
Ranges  which  either  possessed  known  non-linearities  in  their  force 
systems  or  flew  at  large  angles  of  yaw  and  it  was  found  that  their 
motion  could  be  very  well  fitted  by  functions  which  were  solutions 
of  the  linearized  equations.  This  seemed  to  imply  that  the  parameters 
of  these  linear  equations  should  be  "average  values"  of  the 
coefficients  of  the  parent  non-linear  equation.  It  is  the  purpose 
of  this  report'  first  to  derive  the  equations  which  relate  these 
average  values  to  the  non-linear  force  terms  and  certain  characteristics 
of  the  motion  and  then  to  demonstrate  the  great  value  of  these 
relations  by  applying  them  to  a  number  of  programs  which  have  been 
fired  at  BRL.  The  success  of  this  technique  more  than  doubles  the 
values  of  free  flight  ranges  for  both  the  ballistician  and  the 
aerodynamicist «  .  . 

Finally  the  extension  of  this  method  to  the  even  more  important 
problem  of  the  prediction  of  yawing  motion  is  described. 
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2. 


DYNAMIC  CONSIDERATIONS 

In  order  to  investigate  the  yawing  motion  of  missiles  moving  at 
large  angles  of  yaw,  the  exact  equations  of  motion  must  he  derived. 

An  important  feature  of  this  derivation  lies  in  the  proper  selection 
of  the  co-ordinate  system  and  the  dynamic  variables.  This  selection 
should  be  made  so  that  'the  resulting  equations  are  as  simple  as 
possible  and  reasonably  compatible  with  the  basic  assumption  that  their 
solution  may  be  approximated  by  a  solution  of  the  linearized  equations* 

In  other  words  the  non-linear  equation  must  possess  the  same  general 
characteristics  of  the  linearized  equation,  i,e.,  order  of  equations, 
number  of  variables,  symmetry  of  variables,  etc. 

For  arbitrary  rigid  bodies  the  best  coTordina.te  system  is  one 
whose  axes  lie  along  the  body's  principal  axes  of  inertia*  Since  the 
mass  distributions  of  most  missiles  are  rotationally  symmetrical  about 

the  longitudinal  axis  pf  inertia,  all  transverse  moments  of  inertia 

j 

are  assumed  to  be  equal  and  all  axes  perpendicular  to  the  longitudinal 
axis  are  principal  axes  of  inertia.  A  right-handed  Cartesian  co-ordinate 
system  with  numbered  axes  will  be  constrained  so  that  its  1-axls  is 
aligned  with  the  missile' 3  longitudinal  axis  and  pointing  forward. 

D  ,  the  angular  velocity  vector  of  the  co-ordinate  system 
relative  to  an  inertial  system,  will  have  components  (n^,  Dg, 
while  the  missile's  angular  velocity  vector,  <r,  will  have  components 
(o^,  t«2,  to^).  This  definition  of  the  co-ordinate  system,  then,  requires 
that  =  u>0  and  D,  =  aL.  If  D,  and  an  initial  orientation  of  the  2- 
axis  are  specified,  the  co-ordinate  system  would  be  completely  determined. 
In  Appendix  A  the  equations  of  motion  are  derived  for  arbitrary  values 
of  and  we  see  that  these  equations  would  be  greatly  simplified  if 

D-,  =0.  For  this  reascjn  we  will  make  considerable  use  of  a  non-spinning 

i 

co-ordinate  system  for  which  the  2-axis  initially  lies  in  the  horizontal 
plane  pointing  to  the  ljeft  and  whose  axial  spin,  is  identically  zero, 

Turning  now  to  the  question  of  dependent  variables  it  may  be  seen 
that  the  major  contenders  are  the  Eulerian  angles  and  the  direction  ' 
cosines.  The  question  of  compatability  how  appears*  Although  the 
linearized  equations  in  terms  of  Eulerian  angles  are  symmetric  in  the 
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components  of  yaw,  the  corresponding  exact  equations  cannot  be. 
This  is  clear  from  the  definition  of  the  angles.  For  this  reason* 
we  will  express  our  equations  in  terms  of  direction  cosines. 


At  this  point  we  must  precisely  define  our  variables.  The  yaw 
angle  is  defined  to  be  the  angle  between  the  missile's  axis  and  the 
tangent  to  its  trajectory.  If  the  missile's  velocity  vector,  u, 
has  components  (u^,  Ug,  u^)  and  magnitude,  u,  the  sine  of  the  yaw 


,  .T22  +  U/ 

angle  is* —  '"u 


and  a  unit  vector  lying  in  the  plane  of  yaw  has 


direction  cosines 


u.. 


1 


-),  Since  we  will 


“2  +u5 


cOhBlder  only  those  missiles  possessing  trigonal  or  greater  rotational 
symmetry,  it  will  be  convenient  to  represent  quantities  in  the  plane 
normal  to  the  missile's  axis  by  complex  numbers.  With  this  in  mind 
we  define  the  complex  yaw  vector,  X,  to  be  a  vector  which  lies  in  the 
plane  of  yaw  and  whose  magnitude  is  the  sine**of  the  yaw  angle. 


tt2  +  lu? 
u 


(1) 


be  needed  in  this  report  and  will  be 


(2) 


As  can  be  seen  from  an  examination  of  Reference  1  the  Eulerian 
angles  also  introduce  considerable  algebraic  complexity. 

In  Reference  2,  u^  Is  used  as  a  characteristic  velocity  for  forming 

dimensionless  quantities  and,  hence,  the  complex  yaw  of  that  re¬ 
ference  is  proportional  to  the  tangent  of  the  yaw  angle.  Although 
this  use  of  u^  does  simplify  the  center  of  mass  relations,  it  also 

introduces  a  number  of  complications  which  our  choice  of  u  avoids. 


The  cosine  of  the  yaw  angle  will 
denoted  by  l. 


u. 


i  = 


u 


In  a  similar  fashion  the  axial  and  transverse  components  of  the  angular 
velocity  may  he.  separated  and  written  in  complex  form. 


V  - 


oyl 

u 


P  « 


(o)2  .+  ioij)d 


(3) 


u 


where  d  is  diameter. 


Finally  we  will  select  our  independent  variable  to  be  distance  measured 
along  the  trajectory  in  caliberB.  If,  this  variable  is  identified,  by  p, 


P  ra 


at 


(4) 


where  t  is  time. 

3  '  '■  . 

The  generalized  linear  expansion  of  the  aerodynamic  force  and 
moment  assumes  that  the  force  and  moment  are  linear  functions  of  X,  y, 
and  their  derivatives  1 in  a  non-rolling  co-ordinate  system.  For  symmetric . 
missiles  this  assumption  introduces  eighteen  coefficients.  If  we  limit 
ourselves  to  only  those  coefficients  having  a  measurable  effect  on  the 
motion^  this  total  reduces  to  ten  and  we  have. the  following  expansion. 


2  2 


f2  +  iF? 


,2  2 

p  d  u 


] 


i  -  *  p  d“  Kjja  i  (5) 

(  +  ivKp,)X  +  IKgU  “  +  1^1 )  |  (6) 

(7)  '' 

Mg.  +  1M?  -  pd3u2  [(  -vKjp  -  iX^)X  r  K^t  -iK^x'V  tfx)l 

i  •  ”  ■  '  ■  i  ,  -  . 

-|-  ,  r  . 

where  (F^,  Fg,  F^)  are  components  of  the  aerodynamic  force, 

(Mi,  Mg,  MJ)  are  components  of  the  aerodynamic  moment, 


M^  »  - ,  pd3u2  vKa 

■  ,  i 

,3. 


(8) 


p  is  air  density, 

n,a 

A  1 

V  =  —  ,  and 


are  dimensionless  aerodynamic  coefficients. 


For  this  report  we  will  only  consider  non-linearities  caused 

k 

by  the  size  of  the  yaw  angle.  Since  symmetry  requires  that  the 

—  2 

aerodynamic  coefficients  must  be  functions  of  XX  -  5  ,  we  will  Only 
consider  such  a  dependency  on  the  square  of  the  sine  of  the  yaw  angle. 

Two  other  force  coefficients  will  be  used  in  this  report.  These  , 
are  based  on  a  resolution  of  the  aerodynamic  force  in  the  plane  of  yaw 
into  components  along  the  trajectory  and  perpendicular  to  it.  If  F^ 

is  the  component  along  the  trajectory  and  F  the  component  perpendicular 

Jj  ,  / 

to  the  trajectory  and  pointing  toward  the  missile's  nose, 

FD-  -  pd2u2  K,j.  (9) 

Fl  =  -  pd2u2  Kl  S  (10) 

where  6  =Vxx  (5  is.  the  magnitude  of  the  sine  of  the  yaw  angle.) 

In  order  to  obtain  relations  between  these  two  new  coefficients 
and  those  defined  in  Equations  (5  -  6),  we  require  that  the  2- axis 
be  in  the  plane  of  yaw  so  that  X  is  real  and  equal  to  -5. 

Fd-J.F1-8F2  (11) 


FL  *  8  i  F2 


Substituting  Equations  (5),  (6),  (9)>  and  (10)  in  Equations  (ll)  (12) 
and  neglecting  the  Kg  and  terms. 


h  -  ^DA  t  6S  K„  (15) 

'■KL  "  “|KDA  (lif) 

From  Equations  (13)  (14)  the  usual  small  yaw  approximations 
follow: 

kd5kua  (15) 

KL=K»--KM  (16) 

Finally  if  is  eliminated  from  Equations  (13-14)  the  following  useful 
relation  results : 

*L  -  h  '  Kl)  ,  (17) 


In  Appendix  A,  the  exact  equations  of  motion  for  the  aerodynamic 
force  and  moment  as  defined  by  Equations  (5-8)  are  derived  (Eqs. 

(A3,  A10,  A19)i  If  the  effect  of  gravity  and  the  variation  of  spin  are 
neglected  (Jg  =  y  =  D  -  0)  they  may  be  written  in  the  following  form 

in  the  non-spinning  co-ordinate  system  (0  =  0)  * 

ui  -  “■ .jdu  i 1  (18) 

*  "r 

v'  =  0  (19) 

X"  +  |h  -  J  -  ivj  X'  +  M  +  iJL'  -  itfrj  X  =  0  (20) 

where  i 

H  -  i  JL  .  JB  +  ^  -2(JH  -  /  V 

-  ^ 
v  =  |v 

m=V2jm 

j 

T  -  i(JL  .  if2  JT)  | 

A  is  axial  moment  of  inertia 

*  .1  ’ 

B  is  transverse  moment  of  Inertia 

!  t 

k^  is  axial  radius  of j gyration  in  calibers 

k2  is  transverse  radios  of  gyration  in  calibers  f 


J. 

1 


m  Is  mass, 

3.  SMALL  YAWING  MOTION  WITH  NON-LINEAR  MOMENT  AND  NO 

DAMPING 


In  this  section  a  very  simple  example  of  non-linear  yawing  motion 
will  be  considered  in  some  detail.  This  is  the  case  of  a  missile 
flying  at  small  yaw  and  acted  on  by  a  cubic  overturning  or  restoring 
moment.  Later  the  effect  of  other  non-linear  forces  and  moments. as 
well  as  the  geometric  non-linearities  will  be  considered. 

The  basic  feature  of  all  the  theoretical  results  of  this  report 
is  the  assumption  that  over  "small"  sections  of  a  missile's  trajectory 
the  non-linearities  present  in  the  equations  of  motion  do  not  cause 
the  motion  to  be  qualitatively  different  from  motion  based  on  the 
linearized  equations.  For  example,  the  yawing  motion  of  a  symmetric 
missile  acted  on  by  non-linear  forces  and  irioments  should  still  be 
epicyclic  when  a  "small"  portion  of  the  trajectory  is  considered.  The 
parameters  of  the  epicycle,  however,  would  probably  be  related  to  the 
size  of  the  motion.  This  assumption  seems  to  be  reasonable  when  the 
non-linearities  are  themselves  "small".  The  experience  of 
ballisticians  as  based  on  actual  free  flight  tests  have  indicated  that 
relatively  long  sections  of  trajectories  and  large  non-linearities  are 
still  "small"  enough  for  this  assumption. 

In  order  to  illustrate  this  point  the  data  analysis  of  a  non¬ 
linear  drag  force  will  be  outlined.  Since  this  is  the  only  non-linear 
aerodynamic  force  which  up  to  the  present  time  has  been  successfully 
handled  by  ballisticians,  this  digression  will  also  provide  a  good 
background  for  later  derivations. 


The  non-linear  dependence  of  the  drag  force  on  the  magnitude  of 
the  yaw  is  very  well  described  by  the  assumption  that  the  drag  coefficient 
is  a  quadratic  function  of  5  where  6' is  the  sine  of  the  total  yaw  angle. 


(21) 
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In  Reference  5  it  is  shown  that  the  usual  drag  reduction  for  flat 
trajectories*  involves  the  fitting  of  the  time  measurements  to  a 
cubic  in  distance. 


2  ' 

t  *  ao  *  al  P  +  a2  p  +  a3  P' 


(22) 


where  p  is  measured  from  the  middle  of  the  observed  trajectory. 

It  is  further  shoim  there  that  under  certain  reasonable  approximations 
the  following  equation  applies  for  drag  force  of  the  same  form  as 
Equation  (2l).  j 

*  “  *Q  +  ^  P  +  lij  JD0  2“  +  a3  p3  +  ^  JD82  X1  <p> 
where  is  velocity  at  midrange  (p  *  0), 

a^  1b  a  constant  related  to  drag  dependence  on  Mach  number, 


(p)  -  \  \  8  dr  dq,  and 


o  o 


D, 


A,3 


i 


If  the  yaw  is  well  described  by  an  epicycle,  then 
101  102 
^  1  +  Kg  e 


x  =  e 


(24) 


where  K. 


K 


Jo 


-  a.p 

e  .  :  J  j  0^  «  0Jq  +  0j 'p,  and 


KjO*  ^j’  ^jo;*  0jp  8X6  real  constants. 

.2-12  2  i(0x  “  02)  i(0o  “  0p) 

.  .  5  =  XX  j=  K±  +  Kg  +  KlK2  (e  ■  *  +  e  ■  :  1 


)  (25) 


Although  the  expression  in  parentheses  may  be  given  more  simply  by 
2  cos(0^ 02),  thii.3  form  will  be  more  convenient  for  the  non-linear 
yaw  equation.  For!  the  case  of  zero  yaw  drag  coefficient,  we  see,  from 
Equations  (22)  and!  (23),  that 


KD  _  K1X>  ~ 


m 


2a- 


(26) 


i  pd" 


A. trajectory  is:  said  to  be  flat  when  the  component  of  the  gravational 
force  along  thej  trajectory  does  not  change,.  The  yaw  drag  treatment  may 
be  extended  to  jtreat  non-flat  trajectories  if  this  is  necessary. 


When  the  yew-drag  coefficient  has  a  measurable  effect,  it  does  not 
usually  affect  the  quality  of  this  cubic  fit  for  segments  of  trajectories 
as  long  as  10,000  calibers.  It  has  been  found  that  in  thiB  case  the 
same  configuration  flying  at  different  average  yaws  will  provide  good 

2a2 

cubic  fits  but  different  value;;  of  -  .  The  problem  is  then  to  find 

al 

the  quadratic  contribution  of  the  yaw-drag  term  of  Equation  (25). 

Since  all  drag  reductions  contain  a  number  of  periods  of  the 

cosine  term,  the  cubic  fit  of  the  time  history  effectively  neglects 

2 

this  component  of  8  .  With  this  in  mind  we  look  for  the  "best"  quadratic 

approximation  to  I  (p)  over  the  length  of  observed  trajectory,  L. 

1 

2 

To  do  this  we  define  an  average  squared  yaw,  8  ,  so  that  the  integral 


I  (p)  -  6  (§-)  dp  is  a  minimum.*  8  is,  therefore,  the  best 


quadratic  coefficient  from  the  standpoint  of  least  squares.  Differentiating 
the  integral  with  respect  to  8^  and  setting  the  result  equal  to  zero. 


P  Ix(p)  dp 


dp 


Note  that  this  average  squared  yaw  is  different  from  the  raore^ tradition- 


-jr  i\  2  2 

ally  used  mean  squared  yaw  6  (K^  +  Kg  )  dp.  Although  the  average 


squared  yaw  has  better  theoretical  justification,  the  numerical  difference 
is  usually  not  important.- 


(*)' 


p  I 


x(p)  dp 


(27) 


where  I^p) 

0  Q 


(K-j2  +  Kp2)  dr  dq 


=  K, 


10 


-^p 


-  1  +  2af.jP 


(a^'r 


.2 : 


+  k. 
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-  **& 


-  1  + 


aagp 


(aa2) 


The  complete  expansion  for  6  may  be  computed  from  the  power  series 
expansion  of  I^(p).  Since  thj^e  non-zero  terms  are  usually  sufficient, 
the  following  expression  for  6  is  easily  derivable. 

<2  ,  .  -}4 

:TfST)  +3T5I) 


62  -  10 


fl  af  1  KL> 

LK10  L  5l3p  + ;  -TOT)  + 1 

p'j  ]_ 

+  k2o'|_5T5I).  +  TTWJ  +  ~m) 


Returning  to  equation  (26), 

J  2ap  .  ~ 

V  =  -75 -a:  -  +  V  8 

range  pd'  1  6 


] 


(27  ) 


(28) 


/v 

2 


In  most  cases  a  good  approximation  for  8"  can  be  obtained  by  taking 
the  first  term  in  the  expansion  of  Eq.  (27'). 


*W  ^KW'+V(K102,+  K2°2) 


(29) 


The  value  of  this  relation  is  clearly  shown  in  Figure  1  where 

j  /2 

"  is  plotted  i  versus  8  for  a  body  of  revolution  fired  at  yaw 

ang2e§eup  to  50°.  jA  further  check  on  this  technique  lies  in  the 

good  agreement  of  the  wind  tunnel  value  of  IC  p  =  1.85  with  the  free 

|  *0 

flight  value  of  1.82.  (See  Table  IV  on  Page  45. ) 
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Free  flight  tests  have  shown  that  the  situation  for  the  yawing 
motion  is  quite  similar  to  this  yaw  drag  case.  Although  the  non- 
linearities  do  not  measureably  change  the  nature  of  the  epicyclic 
motion  which  is  predicted  by  the  linear  theory,  they  do  affect  the 
values  of  the  epicyclic  parameters.  It  is,  therefore,  necessary 
to  obtain  relations  of  the  same  type  as  Eq.  (29). 

Since  the. simple  case  of  the  small  amplitude  yawing  motion  of  a 
missile  acted  on  by  a  cubic  static  moment  will  first  be  considered. 
Equation  (20)  can  be  considerably  simplified.  The  small  amplitude 
assumption  implies  the  approximations  i  =  1  and  V  =  0  while  the 
restriction  to  a  static  moment  eliminates  the  ferns  in  H,  J'  and  T. 
The  equation  under  consideration,  therefore,  becomes 


62  -  U. 

The  solution  to  the  linearized  form  of  Equation  (30)  for  a  . 
gyroscopically  stable  missile*  is  an  epicycle  without  damping. 


r-~r 


(31) 


■  ■  ■  _  2 

*  A  missile  is  gyroscopically  stable  if  v  >4l^  •  All  statically  stable 

missiles  (Mq<0),  for  example,  are  gyroscopically  stable.  Statically 
unstable  missiles; (^ >^0)  are  gyroscopically  stable  when  their 

_  2 

gyroscopic  stability  factors,  s  =  v  ,  are  greater  than  unity. 
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where  *  :  0  =  0^  +  0^  p 

0-L  02>  311(1 

Ki0'  0^q,  and.  0^  are  real  constants. 

In  Equation  (31)  the  K. n  and  %, n  s  depend  on  the  initial  conditions  hut 

*Wf  ^  ^ 

the  frequencies,  0^  are  functions  fo  the  coefficients  of  the  linearized 
equation  and  are  independent  of  the  initial  conditions.  The  well  known 
relations  for  the  frequencies  are 

'V*  -i  ,  _  v 


\  *%  -  V 


It  will  he  shown,  however,  that  the  frequencies  of  the  quasi-linear 
solution  of  Equation  (30)  do  depend  on  the  initial  conditions. 

If  the  solution  to  Equation  (30)  is  assumed  to  be  of  the  same 
form  as  Equation  (31),  it  can  be  substituted  in  Equation  (30)  to 

provide  the  following  equation 

i0nr  to  -i  p  p  I 

K,„  «  1  -  A  *  V  0,  -  M„  -  M0  (K,„  +  2K.n  ) 


+  K20  e 


r  i  p  t.  p  p  *”1 

3e  +  7  0!  -  M0  -  ^  (K10 

3  b102[-  ^  *  V  0'  -  Mq  -  H,  (K202  t  2K10S)^J  (54) 

f  2  l(2^l  '  0?^  2  i(202  '  M 

|_K10  K20  6  +  K10  E20  6  J  "  °* 


*  The  tilde  superscripts  are  used  to  emphasize  that  the  ^  s  appear 
in  the  solution  to  the  linearized  equation. 
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Unfortunately  Equations  (35  -  36)  are  both  quadratic  equations  and  two 
different  values  of  each  frequency  are  possible.  This  difficulty, 
however*  can  be  resolved.  Since  we  are  considering  solutions  close 
to  the  solution  of  the  linearized  equation,  the  solution  of  Equation 
(35)  close  to^^  should  be  selected  and  similarly  for  Equation  (36).  . 
This  means  that  the  larger  root  of  Equation  (35)  and  the  smaller  root 
of  Equation  (36)  should  be  used. 

p  1 

Equations  (55  -  36)  show  that  Sei,  the  effective  value  of  the 

squared  yaw  for  the  i  -  th  frequency,  is  twice  as  sensitive  to  the 

amplitude  of  other  frequency  as  it  is  to  its  own  amplitude.  The  effect 

comes  from  the  cosine  term  which  was  omitted  from  the  yaw  drag  analysis. 

2 

For  the  yawing  motion  this  periodic  part  of  S  can  not  be  neglected. 

In  order  to  obtain  relations  similar  to  Equations  (32--  33),  M 

_  '  Q 

is  first  eliminated  between  Equations  (35  -  36)  and  then  v  is  eliminated 
between  them 


2  2 

K  -  K 

^+0‘=V  +  M2(^r - — ) 

-0; 

(37) 

f<i  •  02  -  Mo  +  t  4 

(38a) 

2  014-0^1 
e  .  "  01-02 
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(38b) 


-  Kj^2  +  Kjq2  +  K10  '  h  ^0 

|  K-K 

At  first  glance  this  derivation  of  our  "equivalent"  solution 

seems  to  be  surrounded  by  an  atmosphere  of  expediency  and  to  be 

resting  on  empirical  assumptions.  Actually  this  is  not  true.  As  we 

shall  now  show  ouir  technique  has  a  good  theoretical  background  and 

requires  a  minimum  algebraic  work  in  comparison  with  other  methods. 

Finally  it  will  be  shown  that  the  technique  is  not  restricted  to  non- 

2 

linearities  proportional  to  8  but  can  be  easily  extended  to  polynomial 
2 

functions  of  8  . 

'  l  ,  , 

The  derivation  of  Equations  (35  -  36)  may  actually  be  considered 
the  first  step  of  an  iteration  procedure.  It  is  certainly  reasonable 
to  assume  for  small  non-linearities  that  the  first  Btep  of  the  iteration 
has  the  same  form  as  the  linear  solution.  The  error  term, 

"2  [Kio2  e  "  *z)  *  ho  ha  l<202  01>]  - 

introduces  mixed  frequencies  which  are  a  characteristic  of  a  non-linear 
equation.  The  next  step  in  the  iteration  vbuld  be  to  assume  a  solution 
of  the  form 

i  10,  10p  1(201  -  02)  1(202  -  0i) 

X  =  K10  e  +  K20  6  +  IC30e  +K40  e 

Substitution  of  Equation  (39)  in  Equation  (30 )  would  provide  four 

t 

complicated  equations  in  terms  of  the  four.  s  and  two  frequencies 
0^  s,  and  our  error  term  would  then  be  of  the  form 

1(301-2 0o)  1(300-20,)  1(40.-30,)  1(44 

h  <h  e  i  *h  e  +  V  +  \« 

•  i  .  '  , 

where  the  A^  are  fifth  order  combinations  of  the  s .  This  process 

may  then  be  further  iterated  to  yield  a  series  expansion  of  the  almost 
periodic  solution  Of  Equation  (30).  Fortunately  as  we  shall  see,  the 
experimental  results  of  this  report  show  that  only  the  first  step  of 
this  process  is  neieded. 


Another  method  of  treating  Equation  (30)  which  may  seem  to  be  more 
elegant  than  the  method  of  this  report  is  that  of  Kryloff  and  Bogoliubpff . * 
It  will  be  shown  that  this  method  provides  the  same  results  as  our  direct 
substitution  approach  and  requires  more  algebraic  work.  Kryloff  and 
Bogoliuboff  move  the  non-linear  terms  to  the  right-hand  side  of  the 
equation  and  place  the  solution  of  the  linearized  equation  in  these 
terms  so  that  they  are  functions  of  the  independent  variable  p.  The 
method  of  variation  of  parameters  is  then  used  to  solve  the  resulting 
inhomogeneous  linear  equation.  In  order  to  solve  the  differential 
equations  for  the  parametric  functions,  the  terms  arising  from  the 
inhomogeneous  term  of  the  original  equation  are  averaged  over  the  two 
periods  of  the  motion. 


As  an  illustration  Of  this  method  we  rewrite  Equation  (30)  as 


X"  -  iv\  -  MqX  =  Mg  52  X  =  Mg  f(p) 
where  f (p)  *  (K1Q2  +  ^q2)^  e  1  +  (K^2  +  2K1Q2)K20  e  lj52 


m 


2  -  0p)  P 

+  K10  K20  @  +  K20  K10  6 


1(2^  -  $£) 


The  parameters  of  this  solution  are  the  magnitudes  and  the  phase 
angles  Differentiating  Equation  (31)  we  have 


\  =  K10  e  +  K2q  e  +  (Kio  +  i ^10)e 


* 


^1 


+  (*20  +  i^2DK20^6 


If  the  last  two  terms  are  set  equal  to  zero, 


(42) 


*  In  Reference  7  i'k  is  shown  that  this  method  is  fundamentally  the  same 
as  that  of  Van  der  Pol .  They  differ  in  the  fact  that  Kryloff  and 
Bogoliuboff  express  the  linear  solution  in  polar  co-ordinate  form 
and  Van  der  Pol  expresses  it  in  Cartesian  form. 
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then  Equation  (4l)  may  be  differentiated  again  to  yield 
''Vi  9  ^^1  rV  I  9  /v»  I  VI 


i"  =  -  01  k2 


-7b  8  * 


+  iTx  (<0  +  iT10  Ve 


1 %  (4o +  1-?2o  W6 1 2- 


Substituting  Equations  (51*  4l,  45)  for  X,  X  ,  X"  in  Equation  (40) 

t  *f\ 

and  solving  for  +  10-^q  K^q  by  use  of  Equation  (42), 


K10  +  i^10  K10 


M2  f(p)  e 

■%) 


-it 


If  the  right  side  of  Equation  (44)  is  now  averaged  over  a  period  of 
^2  “  ^i*  the  result  divided  into  real  and  imaginary  parts, 


ho  r  0  r 

rt  ..  Kio2  +2K2o2 
?1° ‘  *  %-% 


By  symmetry, 


K20- 


K20  +  2K10 


ho-  “2  -T7fr 


L  1  2  J 

Equations  (45  -  47)  show  that  the  first  approximation;  to  the 
solution  of  Equation  (50)  is  an  epicycle  with  no  damping  but  with 
frequencies  which  differ  from  those  of  the  solution  of  the  linearized 
equation.  The  corrections  to  the  frequencies  are  the  7i0  s  given  by 
Equations  (46)  and  .  (48). 


K102  +  2K202 


0’  =T  -  M  32-——— 
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In  order  to  compare  these  results  with  Equations  (37  _  38) ,  Equation 
(48)  and  (49)  are  first  added  and  then  multiplied.  After  making 
use  of  Equations  (32)  and  (33),  the  following  relations  may  he 


obtained. 


* 


+02 

£•4 


2  2 

t-  .  »  ,10  "  20  \ 

V  +  ) 


r 


%  +M2 


2  2  K10  ^1 
K10  +K20  + 


(50) 

***  K  **»' 
*20  *°2 


X-TC 


-  (Mg  Y 


(51)-. 


(Kio  +  2K2q  )(Kg0  +  2K1q  ) 


20  ,v 20 

^  -%)2 


But  if  second  order  quantities  are  neglected,  Equations  (50  -  51)  are 
the  same  as  Equation®  (37  -  38).  ThuB  for  twice  the  work  we  get  the 
same  result. 


As  a  further  basis  for  our  quasi-linearization  technique  a  special 

form  of  Equation  (30)  will  be  considered  for  which  an  exact  solution 

1b  known.  This  case  is  that  of  a  statically  stable  non-spinning  missile 

16 

in  planar  yawing  motion.  If  X  is  replaced  by  6  e  where  6  is  the 
orientation  angle  of  the  plane  of  yaw,  Equation  (30)  assumes  the  form 


6"  -  (MQ  +  Mg  52)6  =  0 


(52) 


where  M^<  0.  For  planar  yawing  motion  5  must  go  through  zero  and,  hence, 
the  two  amplitudes  must  be  equal  (K^q  =  KgQ  =  K).  Finally  according  to 
Equation  (37)  the  frequencies  for  a  non-spinning  missile  differ  only  in 

ff. 

sign,  (0^  =  -02  =0).  The  quasi-linear  solution  of  Equation  (52)  is 


a  special  form  of  an  epicycle, 


10x  10£ 


S  =  K  (e  +  e  )e 
=  Sm  cos  (0  P  +  -0Q ) , 


-ie 


(53) 
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! 


for  which  Equation  (38)  reduces  to 

where  8m  '«*■  2K.  (maximum  value  of  8) 


(54) 


10 


-©!«=©-  0, 


20* 


The  exact  solution  of  Equation  (52)  is  an  elliptic  function  with 
period  determined  by  complete  elliptic  integral  of  the  first  type.  Thus 
by  use  of  a  little  algebra  and  a  table  of  complete  elliptic  integrals 
it  is  possible  to  compare  the  periods  predicted  by  Equation  (5^  with 
the  period  of  the  exact  solution  of  Equation  (52).  In  doing  this  two  cases 
must  be  considered; 

1.  A  moment  which  grows  faster  than  a  linear  moment  (Mp<0) 

2,  A  moment  which  grows  slower  than  a  linear  moment  (M£?0)  and, 

v  K  5  I 

therefore,  actually  changes  sign  for ■  |  — ^ - 1  =  1. 


Although  a  quasi- linear  approximation  may  be  reasonably  good  for 

,2 


reasonably  large  values  of 


cannot  be  good  for  values  of 


Mg  6 

“2 


«0 


when  Mg  is  negative,  it  certainly 


near  unity  when  Mg  is  positive. 


With  this  in  mind  we  may  now  state  the  surprising  results  of  this 
comparison  with  the  exact  theory.  For  negative  Mg's,  Equation  (82) 
predicts  the  period  with  less  than  1$  error  whenever  the  non-linear  c2 


I  ”2 


8 


M. 


moment  contributiorij  is  less  than  five  times,  the  linear  moment  ( 

.  I  **o 

For  positive  the  error  will  he  less  than  a  percent  when  the  non¬ 

linear  moment  contribution  is  less  than  one-half  the  linear  moment 
2 

I  i 

In  this  case  the  error,  however,  rises  quite  rapidly 


<5). 


( 


Vi 

M„ 


8 


Ki  >■ 


for  larger  angles,  i  As  a  result  Of  these  facts  it  is  reasonable  to  make 
use  of  Equations  (35  -  38)  with  considerable  optimism. 
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Turning  now  to  the  question  of  more  general  non-linearities  an 

examination  of  the  algebra  used  in  Equations  (35  -  38)  indicates  that 

2 

any  polynomial  function  of  6  could  be  used.  In  order  to  obtain 

2 

the  effective  values  of  higher  powers  of  5  we  look  for  those  terms 

i0i 

in  52n  which  upon  multiplication  by  X  yield  terms  in  K1Q  e  or 

i02 

K2q  e  .  Clearly  the  only  terms  which  have  this  property  are 
constants,  e^2  "  $l\  and  e^l  ~  $2) .  In  the  binomial  expansion 
of  &2n  the  only  term  in  which  these  exponentials  appear  are  various 

i(0p  -  02)  ^ 2  "  ft.) 

multiples  of  the  cosine,  i.e.  e  +  e 

Although  the  constant  terms  are  unaltered  in  all  three  forms  of  effective 
yaws,  [&2n]ei,  PL  8114  : [b21^]© ,  the  coeine  term  Joa-kes  different  contri¬ 
butions  to  these  forms.  For  I  &2nJ  ,  its  contribution  }.s  CK?o^ >  for  J~6  |  , 

*  2  r2ri  0X  -  K20d  02  L  -1 

it  is  CK,n  ,  and  for  [6  J  it  is  C  ( - ; - -5 - )• 

01  -  02 

As  an  example  of  this  algorithm  we  calculate  the  three  effective 

4 

forms  of  o  . 

From  Equation  (25),  .  2 

.4  r  .  - hK 


=  a  +  b(e 


2  2 

where  a  =  K1Q  + 


b  =  K10K20 


The  exponential  form  of  the  cosine  allows  an  easy  selection  of  the 
constant  and  cosine  terms  of  the  expansion.  These  particular  terms  of 

4 

the  expansion  will  be  identified  by  brackets  on  5 

[>1  =  a2  +  2b2  +  2a  I  1  (56) 


where 


=  b(i 


=  a2  +  2b2  1  2a 
%  -  02)  1(02  -  0X) 


r  ii  4"]  H4-  2  2 

In  order  to  obtain  [_&  JelJ  °r  l^Je^  K20  *  K10  ’  °r 


K10  ^1  "  K2Q  ^2 

<  -  i 

Equation  (56). 


respectively  should  be  inserted  in  the  empty  brackets  of 


2  2 

=  ,  a  +  2b  +  2a 
el 

_K202] 

2  2 

0  =  a  +  2b  +  2a 
e2 

M 

:  2  2 
=  a  +  2b  +  2a 
e 

X10  ' 

01  l2.*  a.8  ♦  2a  -1  <»> 

L  1*1-4  J 

2n 

In  Table  I  values  of  S  are  given  for  all  values  of  n  between 
one  and  eight. 

i 

4.. .  GEOMETRIC  NON-LINEARITIES  OF  LARGE  YAWING  MOTION 

The  basic  quantities  measured  in  free  flight  are  the  co-ordinates 
of  missile's  center  of  mass  and  the  direction  cosines  of  its  axis  of 

5 

symmetry.  Since  the  usual  formulas  relating  the  direction  cosines 
and  the  motion  of  the  c.m.  to  the  complex  yaw  \  were  derived  for  small 
yaws,  the  exact  relations  have  to  be  derived  for  this  report.  In  this 
derivation  we  will  find  it  convenient  to  keep  the  2-axis  in  the  horizontal 
plane  (fixed-plane  co-ordinates).  Unfortunately  these  co-ordinates  are 
not  the  same  as  the  non-rolling  co-ordinates  (v  =  0)  and  it  is,  therefore, 
necessary  to  calculate  v1  for  these  co-ordinates.  Surprisingly  enough  in 
Reference  8  it  is  shown  that  not  only  is  y  finite  but  it  has  a  non- zero 

'  I 

average  value.  In  our  development  the  cumbersome  Euler i an  angles  of  that 
report  will  not  be  used  and  the  desired  results  will  be  obtained  in  a 
somewhat  simpler  fashion. 

Co-ordinates  in  the  free  flight  range  system  will  be  identified  by 
(xp  Xg,  x^)  and  in  the  fixed-plane  system  by  (y^  y2,  y The  range 

system  has  a  1-axis  pointed  downrange  along  the  intersection  of 

i 

horizontal  plane  and  yertical  plane  containing  the  gun.  The  2-axis  lies 

i 

in  the  horizontal  plane  pointing  to  the  left  and  the  5-axis  up.  The 

fixed-plane  co-ordinates  have  the  1-axis  along  the  missile's  axis,  the 

I 

2-axia  in  the  horizontal  plane  pointing  to  the  right,  and  the  5-axis 

down.  Finally  the  non- spinning  co-ordinates,  which  are  our  fundamental 

co-ordinates  in  the  theory,  have  the  1-axis  along  the  missile's  axis, 

the  2- axis  initially  pointing  to  the  right  in  the  horizontal  plane 

A  ! 

hut  moving  so  that  v  is  zero  and  the  5-axis  fixed  by  the  right-hand 
rule.  ! 
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In  order  to  calculate  the  complex  yaw  it  is  necessary  to  know  the 
components  in  the  range  system  of  the  unit  vectors,  e^,  e^,  along  the 

fixed- plane  axes.  Since  the  fixed  plane  1-axis  lies  along  the  missile's 
axis  of  symmetry,  components  of  the  unit  vector  are  (n  ,  n~,  n, )  where 

:  l  d  j 

the  n^  are  direction  cosines  of  the  missile's  axis  with  respect  to  the 
range  system.  The  restriction  of  the  fixed  plane  2-axis  to  the  horizontal 
plane  is  equivalent  to  the  requirement  that  the  third  component  of  ^ 
be  zero.  This  together  with  the  requirement  that  it  be  perpendicular 
to  the  1-axis  and  pointing  to  the  right  completely  determine  the 
components  of  the  unit  vector  along  the  2-axis. 


According  to  the  right  hand  rule,  the  third  unit  vector  is  equal  to  the 
cross  product  of  the  first  two. 


With  this  information  the  matrix  equation  for  changing  from  the  range 
coordinates  to  fixed- plane  co-ordinates  can  now  be  written: 


30  ■ 


If  the  vector  (4^,  x^,  x^)  is  inserted  on  the  right  side  of  Equation  (63) > 
the  left  side  would  yield  the  components  of  the  velocity  vector  in  the 
fixed- plane  co-ordinates.  These  co-ordinates  will  be  denoted  by  the 
Symbols  (u^,  and  the  symbols  without  the  circumflex  will  be 

reserved. for  the  non- rolling  co-ordinate  system. 


3)  T 


Dividing  by  the  magnitude  of  the  velocity  vector, 

2l  1/2  ' 


r  //ax2\ 

usii  L1  A%;/ 


pc"/  _T  and  simplifying,  we  can  obtain  the  following  expressions 
for  the  components  of  the  complex  yaw  in  the  fixed- plane  system: 


For  the  case  of  a  flat  trajectory,  the  derivatives  are  the  same 
order  of  magnitude  as  n^  and  n^  or  smaller  and,  if  we  make  use  of  the 

o  2  2 

fact  that  =  1  -  -  n^  ,  we  can  obtain  the  first  order  approximation 

of  the  yaw  components 
-  .  A  <2*2 

:X2  "  n2  “  dxx 


(68) 
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Turning  now  to  the  problem  of  computing  the  axial  angular,  velocity 
of  our  fixed-plane  co-ordinates,  we  first  must  state  the  definition  of 
the  angular  velocity  vector  in  terms  of  unit  vectors  along  these  co¬ 
ordinate  axes  e^. 


where  (  )'  denotes  ~  . 


Substituting  in  Equation  (70)  the  co-ordinates  of  the  vectors 

and  from  Equations  (61-62)  we  have 

i 


Similarly, 

A 

M  = 
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A  l 
-  ^  ( 


1  -  n. 


If  the  angle  between  the  missile's  axis  and  a  vertical  axis  is 
denoted  by  ty,  then  n^  =  cos  t  and  Equation  (74)  reduces  to 

v  =  -  fi,  cot  ir  (75) 

Since  is  the  angular  velocity  of  the  missile's  axis  about  the  2-axis 
which  is  fixed  in  the  horizontal  plane,  its  integral  should  be  related 
to  \|f,  Integrating  as  given  by  Equation  (73)>  ve  see  that  this  is  true, 


P2  dp  =  arccos  n^(p)  -  arccos  n^(o) 


t  =  A  -  dp  + 


Although  Equations  (75) ‘and  (77)  are  very  useful  for  the  numerical 
integration  of  the  complete  equations  of  motion,  we  will  find  Equation 
(72)  to  be  much  more  amvenient  for  the  purposes  of  thiB  report. 

If  the  flat  trajectory  approximation  x^  =  pd  is  used,  a  first  order 
approximation  of  v  may  now  be  obtained  from  Equations  (67),  (68),  and 
(72).  ■  , 

C  =  n3„'  -  d3  +  3)  (K  *  -|)  • 

A  At  (78) 

=  X  X  +  T 

»  3  2  , 

y  !  Y  11  X  X  M 

where  T  =  —  X  +  -5-  ^  +  —3—  • 

d 


A  good  approximation  to  the  yaw  components  is  an  epicycle. 
A  A  A. 


"3  10 


10 

cos  $ 

+  K20  cos  02 

(79) 

A 

A 

(80) 

10 

sin  0^ 

+  K20  sin  02 
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where  K^q  are  constants  and 


f  i 

0^  are  linear  functions*  of  p. 

i 

■ '  •  *5*2  =  (K10  hin  +  K2q  sin  $g)(-  K10  ^  sin  ^  -  KgQ  sin  \)! 

2  A.  i 


1 

2 


K10  ^1  +  K20  ^2 


K1q2  cos  2^  +  K202  $gabB2^ 


-  (Ki  +  ^K10K20  sin  sin 

A 


(81) 


of  -  | 


2  ^ 

f  „  2  6* 

Sk 

o 

1  +  K20  ^2 

From  Equation  (8l)  it  can  be  seen  that  v  has  a  non- zero  average  value 

(Fortunately,  for  small  yaws,  this  can 

be  neglected, ) 

Turning  to  the  small'  trajectory  term  T  we  consider  only  the  effect  of 
gravity  and  lift  force  and  obtain  the  following  relatione  for  constant 
velocity  u^  from  Reference  5 ■ 


r  =  Qi  +  JL[jp  sin  \ +  f2-  Bln  %\ 


T  =  ^2  -  ~2  P  -  JL 


K10  ft  .  K20  ft 

__  cos  0  +  -J-  cos  0 

^  % 


LK 


(82) 


(83) 


where  Q  are  small  constants**,  and 

| 

g  is  the  acceleration  due  to  gravity 

pa5 

JL  =  V  KL 


A  little  algebraic  manipulation  results  in  the  relation  that 


* 


** 


The  circumflex  on  the  ^ ’ s  is  to  indicate  that  these  angles  are 

measured  in  the  fixed  plane  system  and  not  in  the  non-rolling 
system  of  Section  2. 

This  is  due  to  the  assumption  of  a  flat  trajectory. 
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T=-(|)  JL2 


K10  iv20 


0 


1 


0 


i 

2  -l 


+  fluctuating  terms  (84) 


*L  6 

Since  for  most  missiles  r,  ^  10”  ,  T  is  effectively  the  same  order 

h 


as  and  we  need  only  consider  the  fourth  order  terms  which  are 

independent  of  the  trajectory  to  obtain  a  complete  fourth  order 
approximation  to  the  non-fluctuating  or  "D.  C.  component"  of  v.  This 
expression  is  calculated  in  Appendix  B.*  If  we  add  T  to  this  and  define 
rtn  to  be  the  "D.  C.  component"  of  v,  we  see  that 


With  these  results  we  can  now  calculate  either  the  exact  position 
or  the  average  position  of  2  and  J-axeB  of  the  non-rolling  co-ordinates. 

Since  our  fixed  plane  co-ordinates  are  turning  at  an  angular  velocity  of 
O'  with  respect  to  the  non-rolling  co-ordinates,  we  have  that 
P 

ij  v  dp 

1  -  X  e  0  (86) 

Ab  it  shall  be  shown  later,  v,  in  Equation  (86)  may  be  replaced  by  |  v J  for 
the  data  reduction  of  large  yaw  firings. 

The  effect  of  large  yawing  motions  on  the  results  of  Section  2  manifests  it 
self  in  two  ways.  First  the  distinction  between  the  non-spinning  co-ordinates 
and  the  fixed-plane  co-ordinates  becomes  important.  According  to  Equations 
(85-86)  these  co-ordinate  systems  rotate  with  respect  to  each  other  with  a  non-z 


*  In  Appendix  B  the  frequency  for  pure  precession  is  compared  with  the 
approximate  relation  based  on  this  expression  for  j"0"|  .  As  can  be  seen 
in  Fig.  15 ,  the  agreement  is  good  for  6  =  sin  45°. 
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averagejjvJ  .  if  0^  and  02  denote  the  epicycle  frequencies  for  the 

range  fixed-plane  co-ordinateB,  it  can  be  seen  from  Equations  (85)  and 
(86)  and  the  assumption  that  the  motion  in  fixed  plane  coordinates  is 
also  epicyclic  that* 


The  second  effect  which  is  the  appearance  of  the  cosine  of  the  yaw  angle, 
may  be  seen  in  the  equation:  of  undamped  yawing  motion  in  the  non- rolling 
system.  This  may  be  obtained  from  Equation  (20)  as 

o  f  _  t 

X"  -  (j  +  iv)x  -  MX  =  0  (89) 


where  M  =  t  k2'2  JM  =  MQ  +  82 


M0  =  k2  V 


*^2  ~  M2  -  2  M0  k2  ^ JMg2  "  2  JMQ^ 


&  .  » 


Equation  (89)  differs  from  Equation  (30 )  by  the  -g  X  term  and  the  presence 
of  the  cosine  term  in  M-  In  order  to  obtain  the^ correct  form  of  Equations 
(37-38)  we  must  derive  a  good  approximation  to  • 
t  /  2  x 1  I 

4  ■  .  (±)  ill  ! 

1  K2J  j2  \ 

!  (^V 

■  -  (5)  — i— 5 

(1  -  »  )  (90) 


For  rapidly  spinning, models  at  high  yaw  Equation  (88).  can  impose  a 

1  J-  -I  _  _  rk  tjj  v)  /W  ^  I,  _ -1  tj-  r^~r  si  f  ** 


large  correction  on  If  an<i  0£  ^ecoiEes  negative 

and  the  product  of  the  rates  in  the  non-rolling  system  would  be  negative. 


36 


=  -  <§)  (B2) 

Inserting  Equation  (25)  in  Equation  (90)  we  have 


1  “5  K10  K20  ^1  "  ^  e 


i(0,  -  0O) 

/  ri*  \  V71  r2'  /j/ 

i  -  0o)  e  +  (Po 


0')  e 


(91) 

l(0o  -  0n  ) 


Substituting  Equation  (91)  in  Equation  (89),  and  operating  on  Equation 
(89)  in  the  same  way  as  we  treated  Equation  (30),  the  f flowing  two^ 
equations  may  be  obtained  from  the  coefficients  of  K^e  1  and  KgQe  ^ 
respectively  0» 

Vj^  +  M0  +  M.*  6^  +  K2O9(0’  -  02)(-|)  -  0  (92) 


22-  +  “0  +  8eS  + 


Eliminating  first  Mq  and  then  v  between  Equations  (92)  and  (93), 


,•  1  2j'  2j' 

(•l  +  0p  +  3  K10  *  *20  *2 


=  v 


*  /K10  ■  K20  v 

7 - 7 - 5 

01  -  02 


j»  ,7  1  2j‘2  2j'; 

■  *1  •  f*2  +  5  K10  <*1  +  K20  *2 


2j'2  *  2 

*2.  Mo  +  “2  6e 


Thus  the  geometric  non-linearity  of  Equation  (89)  introduce  correction 

terms  to  the  left  sides  of  Equations  (37  -  38)  of  the  previous  section 

*  1 

and  replaces  Mg  by  Mg  «  g  Mq  in  those  equations .  If  the  frequencies 

in  Equations  (94  -  95)  are  replaced  by  their  fixed  plane  values  from 
Equation  (87  -  88),  and  fourth  powers  of  K^' s  are  neglected.  Equations 
(94  -  95)  may  be  written  in  the  following  useful  form. 

k*k-s  Lx + oa  - » ♦  *  (Kvy  )  (96) 

L  -1  .  -  02 


5  <Kio2  +  k202> 


*  2 

Mo  +  «2  6e 
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An  examination  of  Equation  (97)  shows  that  the  effect  of  a  large 
percentage  change  in  the  frequencies  implied  by  Equations  (87  -  88)  is 
almost  completely  cancelled  when  the  terms  in  the  cosine  of  yaw  are 
properly  handled.  As  a  check  on  the  derivation  of  Equations  (96  -  97)  they 
are  derived  in  Appendix  C  directly  from  the  yaw  equation  expressed 
in  the  fixed-plane  system.  This  derivation  shows  that  the  assumption 
implied  in  Equations  (87  -  88)  is  equivalent  to  the  basic  assumption 
underlying  this  whole  report. 

5.  EXPERIMENTAL  RESULTS 

In  Section  3;  a  very  simple  example  of  non-linear  yawing  motion  with 

two  degrees  of  freedom  is  considered  in  great  detail.  The  very  same 

i 

quasi-linear  assumption  may  be  applied  to  a  much  more  complicated  system. 

In  appendices  D  and  E  the  yawing  and  swerving  motion  of  a  missile  for 
which  all  aerodynamio  coefficients  are  quadratic  functions  of  5  is 
considered.  Equations  relating  the  parameters  of  the  epicycle 
approximation  (frequencies  and  damping  exponents)  to  the  coefficients 
of  the  parent  non-linear  equation  and  the  amplitudes  of  both  modes  are 
derived.  Since  these  relations  are  essentially  approximations,  their 
value  must  be  determined  by  either  numerical  or  experimental,  checks. 

A  good  check  on  | these  approximations  can  be  obtained  by  means  of 

the  actual  free  flight  motion  of  a  missile  acted  on  by  non-linear  forces 

and  moments.  For  this  reason  various  results  of  firing  tests  on  BRL* s 
9  : 

spark  ranges  were  examined  for  possible  verification  of  the  quasi-linear 
hypothesis.  The  results  of  this  investigation  were  quite  encouraging. 


The  effect  of  a  cubi c  static  moment  on  three  different  programs  was 

first  considered.  These  programs  were  a  finned  missile  program  (8lmm 
10  ! 

mortars  ),  a  body  of  revolution  program  (Army- Navy  Spinner  Rocket 

program^),  and  a  large  yaw  body  of  revolution  program  fired  by  E. 

15) 

Roecker.  ' 

In  the  finned  missile  program  a  free  flight  range  value  of  was 
computed  from  the  prokuct  of  the  frequencies.  According  to  Equation 
(38a)  and  the  definitions  of  M. 


range 


2 =  +  v  6€ 

0  5 
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Thus  the  range  value  of  K^,  as  obtained  from  an  epicycle  fit  of 
individual  rounds,  should  be  a  linear  function  of  the  corresponding 
effective  squared  yaws.  In  Figure  2  the  range  values  of  JC,  are 

p 

plotted  against  and  different  lines  for  each  configuration  tested 
are  drawn  through  corresponding  points.*  The  linearity  exhibited  by 
these  experimental  points  constitutes  the  first  check  of  the  theory. 


It  should  be  emphasized  that  in  all  cases,  epicycles  with  damping 

were  fitted  to  the  data  and  although  the  exponentially  varying  K, ' s 

2 

actually  should  have  appeared  in  the  calculation  of  6g  they  were 
approximated  by  their  midrange  values  K^q's.  The  approximation  is 
implied  by  our  basic  quasi-linear  assumption. 


The  Army-Navy  Spinner  Rocket  program  consisted  of  three  model 

lengths  each  with  three  different  center  of  mass  positions  and  fired 

at  three  Mach  numbers.  A  maximum  of  twenty- seven  values  of  2  was 

& 


therefore  possible.  Since  for  each  case  at  least  four  points  with  a 

p 

reasonable  qpresd.  in  are  .necessary  for  a  good  determination  of  K^2, 

it  was  actually  possible  to  obtain  only  sixteen  values  from  a  careful 
analysis  of  the  126  record  rounds  of  the  program. 

These  values  with  their  standard  errors  are  listed  in  Table  II . 

A  sample  plot  of  the  experimental  points  for  the  9-caliber  model  at 
Mach  number  of  1.8  is  given  In  Figure  J, 

For  some  of  the  Mach  numbers  and  configurations,  wind  tunnel  data 
taken  by  R .  Krieger  were  available .  The  data  for  the  overturning 
moment  were  fitted  by  cubice  and  the  cubic  coefficients  (K^2)  are 
listed  in  Table  II.  At  all  seven  points  of  comparison  the  agreement 
is  good. 


Since  these  finned  missiles  were  not  spinning,  07  =  -  0/  and  the 
effective  squared  yaw  assumes  the  concise  form  ^  (K-^q  +  K20  * ^ 

In  Reference  10,  the  K  *s  were  erroneously  plotted  against 

M 

range 

“p  p  p 

&  i  kiq  +  K20  and,  hence,  the  slopes  obtained  in  that  report  are 
actually  |  KJ  2- 
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TABLE  II 


V 


FOR  THE  ARMY-WAVY  SPINNER  ROCKET 


Mach 

Number 

Forward 

Center  of  Mass 

Middle 

Center  of  Mass 

Rear 

Center  of 

Mass 

Range 

Wind 

Range  Wind 

Range 

Wind 

Tunnel 

Tunnel 

Tunnel 

Five 

Caliber  Long  Models 

1.3 

-  2 

-3  +  1  0 

1  .+  2 

3 

1.8 

-14 

-  7 

3  +  1 

0 

Seven  Caliber  Long  Models 

1.3 

-  21  + 

1 

-  .  5  +  1 

7+2 

5 

1.8 

-  12  + 

1 

-ii 

-  3 

0  +  3 

3 

2.5 

-  28  + 

1 

-26 

-  11 

Nine 

Caliber  Long  Models 

1.3 

-  35  + 

1 

.1  +  .01  * 

1.8 

-  3°  + 

4 

-16+1 

-  3-7  +  -3 

2.5 

-  30. +  2 

-  l  +  .6 

* 

The  center  of  maas  for  these  RCM  models  is  .2  cal.  rear  of  its 
usual  location. 


! 

i 

i 

i 
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Since  models  with  three  different  center  of  mass  locations  were 
fired,  a  second  check  is  possible.  If  the  normal  force  is  expanded 
as  a  Cubic  function  of  &,  the  usual  center  of  mass  relations  provide 
that 


V4)  ■  \ 


^2^  ^2  +  4  ^2 


(100) 


where  (q)  are  the  moment  coefficients  for  a  center  of  mass  located  q 
calibers  rear  of  the  center  of  mass  for  the  *s. 


This  means  that  the  and  ^ 


1  s  for  different  center  of  mass 


locations  are  linear  functions  of  location.  It  was  possible  to  measure 

three  values  of  K  for  only  one  configuration  and  one  Mach  number. 

M^2 

In  Figure  4  both  Kj^  and  2  for  this  case,  are  plotted  against  center 

of  mass  location.  The  fact  that  the  cubic  coefficients  as  well  as  the 
linear  coefficients  fall  on  a  straight  line  is  another  point  in  favour 
of  the  theory.  The  slopes  of  these  lines  are  K_r  ?  and  K-,  respectively. 

For  six  other  cases  it  was  possible  to  compute  2  from  two  values 

of  1C.  g  and  all  seven  values  together  with  the  corresponding  linear 
6 

values  (IC.  )  are  tabulated  in  Table  III.  Once  again  Krleger's  wind 
Uq 

tunnel  data  were  analysed  and  it  was  found  that  at  all  points  of 
comparison  the  wind  tunnel  results  were  in  good  agreement  with  center 
of  mass  values. 

In  Appendix  D,  it  is  shown  that  the  lift  force  and  Magnus  force 
coefficients  may  be  directly  measured  from  the  swerving  motion  when  it  i 
large  enough.  According  to  this  appendix,  most  of  the  swerve  is  associ¬ 
ated  with  the  lower  frequency  and,  hence,  the  range  values  of  and  Kp 

2 

for  each  round  should  be  plotted  against  the  corresponding  value  of  &e0 
A  sample  plot  of  this  process  is  given  in  Figure  5* 


TABLE  III 


Mach 

Number 

V 

' 

Center 

Swerve 

Wind  Center 

Swerve 

Wind 

of  Mags 

Tunnel  of  Mass 

Tunnel 

Five 

i 

Caliber  Models 

1.3 

•  98 

i 

1  -99 

.92  7 

5.2 

+  .2 

5 

1.8 

1.13 

;  1.13 

1.10 

10 

+  2 

,  7 

Seven  Caliber  Models 

1.3 

1.02 

•  98 

20 

10 

+  1 

_  _ 

1.8 

1.13 

1.13 

1.08  12 

9 

+  4 

10 

2.5 

1.21 

'1.20 

1.21  18 

27 

1  k 

18 

Nine 

Caliber  Models 

1.3 

1.06 

1.07 

IT 

17 

+  l 

1.8 

1.14 

1.25 

14.2  +  .4 

15 

+  2 

2  r.5 

1.31 

1.20 

34  ” 

30 

+  2 
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In  order  td  calculate  Kjj  ^  from  2>  necessary  to  expand 

Equation  (17)  in  powers  of  B  and  compare  coefficients 


vv\ 


(101) 


V  ■  V  ■  5  \ +  V 


(102) 


UsingJthese  relations  it  was  possible  to  calculate  eight  values  of 

both  IC.  and  p  from  the  swerving  motion  and  these  are  given  in  Table 
w0  ^B^ 

III.  With  the  exception  of  two  values  of  g  for  the  7“calil>er  models 

6 

-the  agreement  is  remarkable.  The  reason  for  these  two  discrepancies 
is  at  present  unknown. 

Although  the  lateral  displacement  due  to  Magnus  force  is  aboqt  one 
tenth  that  due  to  lift,  it  was  possible  to  make  two  measurements  of 
Kj,  2 ■  The  experimental  points  for  these  two  cases  are  given  in  Figure  6* 
Thi  excellent  internal  agreement  of  these  quite  delicate  measurements 
is  extremely  gratifying. 

Finally  in  Roecker*s  large  yaw  program  it  was  possible  to  check  the 


treatment  for  the  geometric  non-linearities.  In  Figure  7 


range 


calculated  from  Equation  (97 ),  is  plotted  against  .  Here  we  see 
that  the  data  is  essentially  bilinear.  Each  line  corresponds  to  a 
cubic  segment  in  the  moment  plane .  If  the  parameters  of  each  cubic 

are  calculated  from  the  slope  and  intercept  of  its  corresponding  line 

2 

in  the  plane,  they  can  be  pieced  together  to  form  a  smooth 

moment  plot,  (Figure  8).  (An  examination  of  the  spark  shadowgraphs 
revealed  the  fact  that  flow  separation  occurs  at  about  21°  and  this 

explains  the  sudden  change  in  the  moment  curve  at  this  point. )  In 

2 

Figure  9 >  is  plotted  versus  8^  and  the  corresponding  lift  force 
plot  Is  given  in  Figure  10. 

Since  wind  tunnel  measurements  for  this  configuration  had  been 
made  by  W.  Buford,  a  comparison  was  possible.  The  data  was  divided 


Into  angles  of  yaw  less  than  21°  and  greater  than  21°  and  pairs  of  cubies 
fitted.  The  resulting  coefficients  are  compared  in  Table  IV.  The 
agreement  is  excellent. 


As  a  result  of  this  rather  spectacular  success  with  the  non-linear 
static  moment }  the  more  difficult  problem  of  Magnus  and  damping  moments 
was  then  considered.  :  In  Appendix  E  it  is  shown  that  a  quadratic  dependence 
of  Kg,  and  Kg  -  on  5  so  affects  the  damping  of  the  epicycle  that  the 
usual  linear  formulas  for  Kg  -  and  Kg,  are  actually  the  following 
combinations  of  coefficients  j 

(kh  *  v 


range 


k  -  % 


range 


;■  “v  p .  p  (Kio2  ■  ^o2) 

j  1  2 

!  t-  2  j  1 2  2ji  2 

V  J  ir*  *  2  .  v  *  r*\  ,ho  h  "  *20  ^2 

*»  V  «  V?)  ( 


) 


(10?) 


(104) 


where 


\z  '  \z  -  ^2  + 1  hm 

^.2  "  ^Te2  '  S;*TL. 


Range  values  of  Kg  -  K^  and  Kg,  for  the  Army-Navy  Spinner  Rocket 

program  were  fitted  by  least  squares  to  Equations  (103-104)  and  values 

0f  ^0'  ^0  “  SlAO'  ^Te2  and  V  2  were  confuted.  In  all  cases  the 
*  6  S 

coefficients  of  g  were  small  and  this  quantity  was  poorly  determined* 
K^2  was  therefore  omitted  from  Equations  (105-104)  and  they  were  fitted 
separately  to  range  vhluee  pf  K  -  and  K^* 

Since  this  measurement  depends  on  the  damping  exponents  and  hence 

1 

is  quite  delicate^  it  ;was  possible  to  obtain  only  eight  values  of  tL,  0 

|  .  *  r52 

from  Kg,  and  four  from  (Kg  -  K^J  .  These  are  listed  in  Table  V 

range  range 
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TABLE  IV 


Large  Yaw  Results  for  M  =  2.3 


V 

0 

V 

h 

Lo 

\  2 

Wind  Tunnel 

1.85 

.825 

-1.1 

1.01 

4.o"l 

before  separation 

Range 

1,82 

.829 

-  1.0 

0.99 

3.6  j 

5  =  sin  21° 

Wind  Tunnel 

.729* 

0.0 

1.45* 

0.0  1 

after  separation 

Range 

.719* 

0.1 

1.43* 

o.4  J 

5  =  sin  21° 

*  Evaluated  at  separation 
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TABLE  V 


Mach 

No._ 

V 

V! 

°V 

V  V 

OJ 

C.M. 

Swerve 

FCM  MCM 

RCM 

Seven  Caliber 

Long  Models 

1.5 

l4 

17  +  1 

4.4 

K^,  -16  +  5-  5  +  1 

_  _ 

'range 

(kh- 

W  -  20  +  14  - 

range 

— 

1.8 

22 

-  - 

4.8 

-34+16 

0  1 

' 

1 range 

<kh- 

V  -  25  i  14  - 

range 

2-5 

19 

5.1 

K  -36+5 

range 

-6  + 

(Kj- 

V.  -26J4  --  ' 

range 

-- 

Nine  Caliber 

Long  Models 

1.8 

1? 

6,5 

K,  -29+2 

range 

-5.2  + 

(Kh- 

K«!  -  26  +  5  - 

range 

— 

*  Kp  2  l°r  Mach  number  1.3  was  given  in  Figure  6  as  8+1. 


NOL  Wind  Tunnel  measurements  were  nja.de  by  Luc  huh  and  Sparks  for  7 
caliber  model  at  M  between  1.6  and  2.5.  According  to  these  results 

V2  “  n,  CP  i  54. 

5  5  i 
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and  the  agreement  between  pairs  of  K_  ?  is  within  the  standard  errors.* 

i5 

In  Figures  11  -  12,  a  pair  of  plots  of  experimental  points  for 

range 

and  (Kg  -  K^)  are  given. 

range  : : 

In  all  four  cases  values  of  2  for  two  center  of  mass  locations 

5 

1  s  were  computed  from  these.  At 

the  Naval  Ordnance  Laboratory,  Luchuk  and  Sparks1"*"  have  made  wind  tunnel 
measurements  of  these  quantities  and  their  results,  as  given  at  the  bottom 
of  Table  V,  are  in  reasonable  agreement  with  flight  tests.  Finally  it 
should  be  noted  in  the  Table  that  at  one  point  of  comparison  of  center  of 
mass  and  swerve  values  of  K_  0,  (  the  seven  caliber  models  at  M  =  1*3)  the 

agreement  is  good. 

As  a  last  example  of  this  quasi-linear  technique  we  will  consider  a 
body  of  revolution  which  displayed  an  extremely  non-linear  Magnus  moment. 

The  range  values  of  K„  for  this  model  were  obtained  by  E.  Roecker  and  are 

2  1 

plotted  versus  5  in  Figure  13.  Although  these  data  are  fitted  by  two 
e 

lines,  only  the  first  line  is  well  determined.  In  fact  if  the  fourth  point 
from  the  right  were  neglected,  the  large  yaw  values  would  be  reasonably 
well  represented  by  the  dashed  horizontal  line. 

In  Figure  lb  the  corresponding  cubic  segments  are  plotted  and  compared 
with  a  BRL  wind  tunnel  curve  obtained  by  A.  Platou.  The  good  qualitative 
agreement  for  this  strongly  non-linear  moment  is  remarkable. 

Since  the  quasi-linear  relations  have  been  so  well  verified  by 
experiment,  they  should  be  quite  valuable  in  the  Important  problem  of  the 
prediction  of  missile  motion  from  a  knowledge  of  the  force  and  moment 
curves.  This  application  of  the  quasi-linear  technique  is  described  in 
more  detail  in  the  next  section. 


were  obtained  and 


v 


'  s  and  CP, 


Thus  an  observed  dependence  of  (K^  —  K^)  on  magnitude  of  yaw  is 

range 

not  necessarily  due  to  non-linear  damping  moments  hut  may  he  due  to  a 
non- linear  Magnus  moment . 
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PREDICTION  OF  NON- LINEAR  YAWING  MOTION 

Although  the  application  of  the  theory  of  this  report  to  range  data 
analysis  enhances  markedly  the  value  of  spark  ranges,  an  even  more  important' 
application  is  the  prediction  of  the  yawing  motion  of  a  missile  acted  on 
by  non-linear  forces  and  moments.  The  experimental  results  of  this  report 
indicate  that  this  inverse  problem  should  be  successfully  handled  by  the 
same  theory.  In  this  section  we  will  outline  the  procedure  for  this 
prediction  problem. 


First  the  sizes  of  the  two  arm's  must  be  obtained  from  the  initial 
conditions 


101  i02 

^0  =K1  6  ,  +K2  6 


(105) 


.  10, 

(-  ax  :+  101)K1  e  x  +  (-  at. 


(106) 


In  order  to  take  care  of  the  effect  of  damping  the  trajectory  should 
be  divided  into  intervals  over  which  neither  amplitude  changes  by  more 
than  50 fyo  The  length  of  the  first  interval  may  be  estimated  from  the 
linear  damping.  If  fhe  calculated  quasi-linear  damping  is  much  larger, 
the  interval  may  then  be  shortened  and  the  process  iterated.  The  lengths 
of  the  other  intervals  should  be  determined  in  a  similar  manner.  Values 
of  0^,  frequencies  in  non- rolling  co-ordinates,  can  now  be  computed  from 
Equations  (92  -  93)  which  are  simple  quadratic  equations.*  (0,/  may  be 
computed  from  Equations  (87  -  88).)  Then  values',  of  niay  be  computed 
from  Equations  (e4  -  E5)  which  are  even  simpler  linear  equations.* 

. 

*  For  ease  of  calculation  the  linearized  val\ies  of  0.s  and  ot^s  may  be  first 
placed  in  the  coupling  terms  arising  from  the  non-linear  geometry  and  an 
iteration  performed  if  needed. 
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If  a  more  complex  polynomial  dependence  of  KM  and  on  &  is 
required,  this  may  be  handled  in  the  following  way. 

Let 


(107) 


K' 


n 


K  5“^ 

^  ' 

.  o 


(108) 


k  =  o 


Then  Equations  (92  -  95)  and  (E4  -  E5)  take  on  the  following  slightly 

* 

more  complex  form: 


n 


4  -  v0; 


d f2 

0o  “  v 


k  =  0 


(&2k)el  +  K20^1  ~  *Zi2)(2£)  +  ai(H  "  0!1)  =  °  (109) 


2 ,  j  ' 


n 


4  +]T4<a2k> 


0 

e2  +  K1Q2^2  "  ^1^2^  +  ^  =  0 


k  =  0 


n 


a 


:p ( 20^  -  y)  -  HO01  +  V 


k  =  0 


4<62k>el  -  V  [<K102  +  *20*  *1  +  K202P2 


-  JL,2^1  ”  ^2^K20  +  ^l^°aK10  +  a2K20  ^ 

o 


(111) 


,  t  ,  t 


+  (- 


;ai^2  +  a2^1w  2 


>K2cT  =  0 


a2(202  -  v)  -  H j'0  +  v 


Or2 


T2k(82k)e2  -  H&2  tlO2  +  “zO^S!  +  “loX 


k  -  0 


Eqs.  (92-93)  are  also  slightly  modified,  by  the  presence  of  damping* 

kQ 


and  those  terms  appearing  in  the  right  hand  summations  which  have  powers 

of  6  greater  than  2 n  are  neglected.* 

! 

It  should  be  remembered  that  these  equations  are  derived  on  the 
assumption  of  a  flat  trajectory.  The  complications  introduced  by  gravity 
Will  be  considered  in  a  later  report.  In  any  event  for  a  large  class  of 
problems  this  simple  procedure  should  allow  reasonably  good  prediction  of 
the  motion  of  missiles  acted  by  non-linear  forces, 

7.  SUMMARY 

1.  A  convenient  expansion  of  the  aerodynamic  force  system  for  large 
yaw  has  been  obtained  and  the  necessary  geometric  relations  for 
these  large  yaws  derived. 

2.  Relations  for  parameters  of  a  linear  equation  which  is  equivalent 
to  the  actual''  non-linear  equation  have  been  derived. 

3.  These  relations  have  been  tested  by  actual  firing  tests  and 

excellent  internal  consistency  has  been  observed.  Where  wind 
tunnel  measurements  were  available,  good  agreement  has  been 
obtained.  j 

ft  i  p 

The  effect  of  the  terms  in  J  which  are  higher  order  than  6  on  the  ' 

t  L 

JL  X.  term  in  Equation  (20)  has  been  neglected. 
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